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Abstract 



, In a first attempt to describe the effect on neutrino travel inside 

the Earth caused by general relativity in the case of a dense Earth, we 
■ have neglected the Earth's rotation, the Earth's ellipticity and also the 

, surface terrain variation, nevertheless we have focused our attention on 

\ the density description of the Earth interior provided by geophysic's 

models such as PREM. Assuming a non rotating Earth, the general 
relativistic effect on neutrino travelling from CERN to GRAN-SASSO 
happened to produce a delay of St = 4. 1863 picosecond. 



X 

£ : 1 Introduction 



After the paper by the OPERA collaboration (cf. pQ) claiming an ad- 
vance of the CNGS neutrino beam from CERN to Gran-Sasso of (57.8 ± 
7.8(stat.) + 8.3 — 5.9(sys.))ns, strong excitement occurred around all pos- 
sible General Relativistic corrections that could account for such a strong 
effect, leaving a big anomaly on neutrino speed faster than the speed of light 
of (v — c)/c — (2.37 ±0.32(stat.)(sys.)) 10~ 5 . Two main corrections are to be 
applied about this measurement: Firstly the Sagnac Effect ([3]) (not detailed 
here), which expresses the impact of the earth's rotation around its axis on 
time synchronization with atomic clocks using GPS in common view mode, 
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and secondly the General Relativistic effect on the neutrino travel due to 
the presence of the mass of the planet Earth. It is this second effect, that 
is deeply detailed here, with the implication of a non homogeneous distribu- 
tion mass density inside the earth. While a simple non-rotating model with 
spherical symmetry with geophysics density data input, a first calculation is 
achieved, providing a 5t = 4.1888 picosecond delay effect. Some second or- 
der effects are further taken into account such as the ellipticity of the earth, 
the effect of mountains and of geoid fluctuations, meanwhile the same delay 
effect of 5t = 4. 1863 picosecond is computed. Of course such an effect is 
far from being measured due to the timing resolution of the most uptodate 
technics of nowadays neutrino detection and also neutrino beam production 
with accelerators timing technics. Previously Kiritsis and Nitti (cf [4]) have 
calculated the General Relativistic effect of neutrino travel between CERN 
and Gran-Sasso, but they assumed a metric of a black hole or far from a 
mass supposed to be ponctual which is not the case here where neutrinos are 
going inside the earth. Of course if this approximation is slighly wrong at 
the surface of the earth and on the crust, it is deeply wrong if the neutrino 
beam is crossing the earth from one side to the other; with how modelization 
of density distribution inside earth, a delay effect of 86. 65pico second is to be 
expected for a full crossing along the earth diameter. The earth is assumed 
not to be rotating, so that a Schwartschild metric is sufficient, while in the- 
ory a detailed calculation in the frame of a Kerr metric would be needed, 
nevertheless, this effect is still extremly tiny. This effect cannot account for 
any advance of the neutrino. Recently the ICARUS experiment also carrying 
measurement on the CNGS neutrino beam did not notice any discrepancy 
between the speed of neutrino and the speed of light (cf |2J). 

2 Toward a relativistic metric inside the earth 

As we do not take into account the rotation of the earth and a spherically 
symmetric density model of the earth interior, we can use the Schwartzchild 
metric to take into account general relativity as follows: 



g = (l- —^)dt 2 - (1 - — — YHt* - r 2 (d9 2 + sin 2 6d<p 2 ) (1) 



The Schwartzchild metric can be easily rewriting using $ potential : 



2GM 



2GM 



g = (l- $)rft 2 - (1 - $) -1 dr 2 - r 2 (d9 2 + sin 2 9d<p 2 ) 



(2) 
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$ is defined by : 



. 2GM 2V 
= — = (3) 



V is easily identified as the classical Newtonian gravitational potential: 



V — VNewton ~ (4) 

r 

At this point, a remark must be made. If Tolman-Oppenheimer-Volkoff equa- 
tion (TOV) is perfectly suited to describe the interior of a spherically sym- 
metric isotropic interior of a neutron star in isostatic equilibrium between 
gravity and pressure with strong general relativistic effects with the input of 
Equation of State (EOS) of the nuclear matter, in the case of the earth the 
situation is slighlty different. The detailed atomic composition of the inner 
earth is a big unknown field and the knowledge of its interior (density and 
pressure) is mainly coming from the long standing experience of geophysicist 
modelling the seismic waves propagation through the interior of the earth. 
Then the TOV equation applied to the earth is directly short cut in this 
work and the density distribution from geophysics is directly plugged into 
the Schwartzchild metric as an input. 



3 Computing the gravitational potential with 
Gauss theorem and Poisson equation 

Because of spherical symmetry hypothesis, the density inside the Earth, 
varies only with r. While differential formulation of Gauss's law for grav- 
ity is given by: 

AnGp = - A V (5) 
The spherical symmetry implies for the expansion of the Laplacian: 

The integral form is given by the Gauss theorem : 

r r 

V(r) -V = -AuG j ^ j p(r)r 2 dr (7) 
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M(r) is defined as the Earth's mass integrated between the center of the 
earth and radius r. 

r 

M(r) = J p(r)Aixr 2 dr (8) 
o 

Gravity acceleration is directly obtained through: 

r 

g(r) = — - / p(r)A-Kr 2 dr 



The classical Newtonian gravity potential is obtained without difficulty as 
follows. 



r 

V(r) = -j 



GM(r)dr 



+ V 



o 

r 



g(r)dr + V (10) 



o 



$( r ) = 



2^(r) 

' 2GM(r)dr 



r 2 c 2 
o 

r 

2g(r)dr 
c 2 

o 



$ (11) 



*o ^ (12) 



The second constant term V and $ m the integral accounts for the prescrip- 
tion of a flat Schwartzchild metric at an infinite distance r from the Earth 
center, as the newtonian gravitational potential is plugged into Schwartzchild 
metric through $, in such a way that $(+oo) = 0. 
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4 The radial density of earth interior obtained 
by geophysics methods with seismic waves 
propagation 

Since many years, geophysicists have been monitoring earthquakes and seis- 
mic waves propagation through the earth interior. Two main kinds of seismic 
waves crossing the Earth are observed, primary waves (faster) and secondary 
waves (slower). Beside the difference of speed of these two categories of waves, 
P-waves are compressional waves, while S-waves are shear waves which do 
not propagate in liquids. The propagation of these waves inside the Earth 
and their sophisticated diffraction processes at the interfaces between solid 
and liquid zone of the inner earth, inform seismologists about the nature of 
the phase of the matter and its density versus depth inside the Earth. These 
two speeds can be calculated as follows: 



^primary 



(13) 



V secondary y P ("^) 

Bulk modulus =K 
Shear modulus p 
Density =p 

A seismic parameter Q S eismic{r) is obtained cancelling the shear modulus out 
of the primary and the secondary waves. 

^ seismic( r ) — ^primary ~ ^secondary (^) 

Moreover, assuming hydrostatic equilibrium and adiabatic exchanges, Adams 
and Williamson (1923) showed a relationship between this geophysic param- 
eter Q S eismic, the density and the gravity acceleration at a distance r from 
the center of the Earth: 

dp = p(r)g(r) , , 

dr $ • ( r ) ^ 1 

^ seismicv ) 

The numerical integration of the Adams and Willimason equation with & S eismic{ r ) 
gives the density (see Figure [I]). Then we computed the gravity acceleration 
versus r (see Figure [2] ) and afterward the Newton gravitational potential 
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inside earth and so the $ potential term of the Schwartzchild metric inside 
the Earth (see Figure [3]). The $ se isrmc( r ) input data are obtained with 
a modelization of the propagation of primary and secondary seismic waves 
for any depthes inside the Earth in order to reproduce the observed seismic 
waves patterns measured by seismographs scattered all over the earth sur- 
face. Each seismograph is recording a seismic wave train characteristic of 
the angle between its position and the epicenter of each earthquake under 
study, and also of the depth of the epicenter below the surface. Nowadays, 
sophisticated models exist (PREM, PEM,IASP91) after the first attempts 
by Bullen (cf. [5]) and others (1937). While the most important part of the 
surface of the Earth is covered by ocean, a continental Model PREMC does 
exist and we have used this model with a density of 2.72tons/m 3 for the 
upper continental crust instead of l.tons/m 3 for the ocean. We focused our 
attention on the Preliminary Earth Reference Model (PREM) by Dziewon- 
ski & Anderson (1981) (cf. [6]) which is highly sufficient for our purpose of 
integrating the density along the radius of the earth in order to obtain firstly 
the acceleration of the gravity, secondly the newtonian gravitational poten- 
tial inside the earth and thirdly a spherically symmetric Schwartzchild-like 
metric inside the earth ( cf. [7], [8j & [9j). (see Figures El [3] H]) 



r = 



R/6371.km 



(17) 



PREM shows various zones inside earth: 



Inner Core (solid): 



O.km <R < \2\7.\km 



0. <r < .19103 
p(r) =13.0122 - 8.45292r 2 
$(r) =2.4845 10~ 9 - 1.6396 10~V 
+ 3.1953 10~ 10 r 4 



(18) 
(19) 
(20) 



(21) 



Outer Core (liquid): 



1217.1fcm <R < 3458.7km 



19103 <r < .54712 
p(r) =12.5842 - 1.69929r - 1.94128r 2 - 7.11215r 3 
$(r) =2.4950 10" 9 - 1.3621 10" 10 r - 1.0558 10~V 



(22) 
(23) 
(24) 



- 1.0388 10" 9 r 3 + 1.3768 10"V - 4.3678 10~ 10 r 5 



(25) 
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Lower Mantle: 



3458.7km <R < 570 Libra (26) 

.54288 <r < .89483 (27) 

p(r) =6.8143 - 1.66273r - 1.18531r 2 (28) 
$(r) =3.3847 10~ 9 - 4.2421 10~ 9 r + 5.1207 l(T 9 r 2 

- 4.0832 10~V + 1.2106 10 _9 r 4 (29) 

Transition Zone: 

5701. km <R < 5951. km (30) 

.89483 <r < .93407 (31) 

p(r) =11.1198 - 7.87054r (32) 
$(r) =2.1071 10~ 9 + 6.0380 10~ 10 r - 1.63960 10 _9 r 2 

+ 3.19532 10" 10 r 4 (33) 

Low Velocity Zone: 

5951. km <R < 6336.km (34) 

.93407 <r < .99450 (35) 

p(r) =7.15855 - 3.85999r (36) 
$(r) =2.7521 10~ 9 - 8.0668 10~ 10 r - 1.4230 10~ 9 r 2 

+ 1.2167 10~ 9 r 3 - 3.4782 10~ 10 r 4 (37) 

Lower Crust: 

6336. km <R < 6351. km (38) 

.99450 <r < .99686 (39) 

p(r) =2.92 (40) 
$(r) =3.5782 10~ 9 - 2.3585 10~ 9 r - 1.4790 10~ 10 r 2 

+ 3.1953 10~ 10 r 4 (41) 

Upper Crust: 

6351. km <R < 6371.km (42) 

.99686 <r < 1. (43) 

p(r) =2.72 (44) 
$(r) =3.6310 10~ 9 - 2.4644 10~ 9 r - 9.4700 10~ n r 2 

- 3.1953 10~ 10 r 4 (45) 
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5 Integrating the travel time of the neutrino 
along a geodesic inside earth 

5.1 Shapiro delay inside earth 

In the sixties, I. Shapiro proposed a test of General Relativity using the delay 
of radar bouncing over Mercury and Venus checking the relativistic effect due 
to the mass of the sun as the radar beam towar these planets was coming 
closed to the sun, this is known as time delay of radar echoes (cf [TU], 
[12j). The Lagrangian of the trajectory of the neutrino inside earth: 

C = l -g^x v (46) 

Inserting the metric into the Lagrangian: 

2£ = (1 - $)i 2 - (1 - <&y l r 2 - r 2 (9 2 + sin 2 9^ 2 )) (47) 

Along the orbit, the Lagrangian C is constant. A simple change of coor- 
dinates can be applied with the condition 9 — ~ which implies : 9 = 

=>- 2C = (1 - $)i 2 - (1 - $) _1 f 2 - r 2 2 (48) 

Applying Noether's Theorem to the invariance of this Lagrangian (and of the 
$ potential) respect to the variables and t , two conserved quantities are 
extracted L and E as follows : 

BC 

— —r = r 2 <f) = L = cte (49) 



dC 

r- = t(l - $) = E = cte (50) 

dt 

After some simple change in the equation, we get: 



(51) 



E 



(52) 
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Neglecting the mass of the neutrino (few milli-eV) respect to its energy (17 
GeV), the Lagrangian is set equal to zero as in the photon case: 

=^2£= (l-$)- 1 (E 2 -r 2 )-^ = (53) 



Simplified as: 

L 2 



r 2 = E 2 - (1-$)— (54) 



r 



The derivation of r respect to the time is obtained here: 



Thus: 



rlr T 2 

Q 2 £ 2 (l-$)- 2 = £ 2 -(l-$)±- (56) 
at r z 

Thanks to the spherical symmetry of this central potential, the neutrino 
crosses the closest point to the center of the earth at a distance r = ro; this 
distance is a minimum, then f = 0. At this minimum distance, the potential 
is also an extremum, this value is assigned to $ = $ . 

^£ 2 -(l-$ )^ = (57) 



« <1> 2 = T^k (58) 

L and E are replaced by r and $o : 

dr , . 

dt = = (60) 

V ; V l-$o r J 

The duration of the travel from the surface (point A) to the minimum dis- 
tance r is calculated through this integral: 

■'■A 

t A -to= [ d \ (61) 



ro 
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The total duration of the travel of the neutrino from point A to point B is 
computed adding the two integrals: 



tB — tA — \tB — to\ + \t A — t \ (62) 
tB — tA 



dr , f dr 



$ r V 1 — $o r 

(63) 

5.2 Numerical application with PREM based metric 

The computation of the general relativistic effect for a neutrino baseline of 
730 km leads to a delay of 4. 1888 picosecond, which is equivalent to a delay 
length of 1.255 mm . In other words the effect decreases the apparent speed 
of light of a relative amount ^ = 1.7199 10~ 9 . This computation is obtained 
with a Monte-Carlo integration method applied to equation (63) with 10 9 it- 
erations. We have computed this effect for 200 distances of neutrino baseline 
crossing the earth spanning between 63 km to 12700 km each one based on 
the same Monte-Carlo method (see figures 0, [6])- The computed result, in 
the hypothetical case that a neutrino beam crosses the whole earth diameter 
(12700km), is a global delay effect of 86.65 picosecond, which is equivalent 
to a length = 2.599cm. This delay length is a one "inch effect". Also the 
apparent speed of light (neutrino speed) is decreased of a relative amount 
^7 = 2.046 10~ 9 . In any case, no advance compared to the speed of light is 
seen here, only a delay through this general relativistic effect due to the gravi- 
tational potential caused by the mass distribution inside the Earth. Anyhow, 
this calculation is done with the following assumptions, no rotation of the 
earth, the earth is supposed to be perfectly spherical and the computation 
is done from one point of the surface to another point of the surface with 
neither mountain, nor any relief, with neither an ellipsoid shape of the earth, 
nor a geoid frame of the earth (i.e. without any angular anisotropic distribu- 
tion of masses inside the earth). A fit of the computed delay 5t[pico second] 
versus the neutrino baseline length b [km] is expressed as follows: 

8t\ps] = .1599 + .55318 10 2 & + .60543 10 7 & 2 + .29928 10~ n 6 3 (64) 
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6 Further improvements beyond this simple spher- 
ical model 



6.1 Effect of the ellipticity of the earth within the WGS84 



The ellipticity of the earth is fully described with the World Geodetic System 
84 (WGS84) (see: [13]), this geodetic system describes an ellipsoid as a good 
approximation of the Earth surface. The origin of the ellipsoid is the center 
of the earth. This ellipsoid is characterised by two primary parameters a and 
e: 

Semi-major axis = a = 63 78137. m 
Eccentricity = e = = 0.081819190842622 

There are also some secondary parameters: 
Semi-minor axis = b = 6356752.3142 m 
Flattening = / = ^± = 0.00335281066474 

Any position on this ellipsoid can be referenced with the longitude, geocentric 
latitude (angle refered to earth center) or geodetic latitude (angle refered to 
the normal of the ellipsoid) and radius, or in cartesian (x,y,z) coordinates. 
Geodetic latitude = <fi 
Geocentric latitude = 0' 

Longitude = A The relationship between the geodetic and the geocentric 
latitude are defined as follows: 



The relationship between the cartesian position (x,y,z) and the latitude, lon- 
gitude plus radius cordinates are given by: 



frame 



tan((f)') = (—) 2 tan((f)) 



(65) 



x 



acos(<f))cos(\) 



(66) 



y 



y/l — e 2 sin 2 ((f)) 
acos(<f))sin(\) 



(67) 



z 



y/l — e 2 sin 2 ((f)) 
a(l — e 2 )sin((f>) 



(68) 



a/1 — e 2 sin 2 ((f>) 




(69) 
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The value of the gravity acceleration at the surface of the ellipsoid is easily 
computed with the geodetic latitude : 




(70) 



Normal gravity at the equator = 7 egu . = 9.7803253359 m/s 2 



7 (0) = (9.7803253359) 



(1 + 0.00193185265241sm 2 ((/>)) 



m/s 2 (71) 



- O.OO669437999O14sm 2 (0)) 



6.2 Positions of the trajectory of neutrino from CERN 
to Gran-Sasso in ECEF WGS84 frame 

In the Earth Centered, Earth Fixed (ECEF) coordinate system, the position 
of CERN reference point ( T-40-S-CERN) of the CNGS neutrino beam is de- 



x = 4394369.327 y = 467747.795 z = 4584236.112 
Latitude =46 U ' 32.793 " N Longitude =6 °4 ' 32.987 " E 
Altitude referred to earth ellipsoid (WGS84)= 429.72 m 
Distance to earth center= 6367455.61 m 

The Opera reference point (Al-9999) is defined as: 

x = 4582167.465 y = 1106521.805 z = 4283602.714 

Latitude = 42 °27' 13.324 " N Longitude = 13 34 1 34.035 " E 

Altitude referred to earth ellipsoid (WGS84)= 1012.22 m 

Distance to earth center= 6369450.50 m 

The middle of the neutrino trip between CERN and OPERA reference points 
is located in the North of Tuscany, in the Parco Nazionale dell'Appenino 
Tosco Emiliano, close to the village of Casalina, near the towns of Potremoli 
and Montelugno : 

Latitude = 44 °24 ' 34.289 " N Longitude = 9° 46' 49.737 " E 
Depth referred to earth ellipsoid (WGS84) = 9740.46 m 
Distance to earth center= 6357969.42 m 
Radius of ellipsoid at the surface (WGS84)= 6367709.88 m 

The point of the closest approach of the CNGS neutrino beam to the center 

of the earth is located near Barzana di Sotto at : 

Latitude = 44 °29 ' 59.289 " N Longitude = 9° 46' 5.733 " E 

Depth referred to earth ellipsoid (WGS84) = 9730.56 m 



fined as: 
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Distance to earth center = 6357945.63 m 

Radius of ellipsoid at the surface (WGS84) = 6367676.20 m 



The coordinates of the position where the CNGS neutrino beam is exiting 
from the earth toward the sky is located on Campo Imperatore not far from 
a crossroad toward Santo Stefano di Sessanio, at a distance of 9750m from 
Opera reference point: 

Latitude = 42 °24 ' 0.278 " N Longitude = 13 40 ' 10.592 " E 
Altitude referred to earth ellipsoid (WGS84) = 1585.49 m 
Distance to earth center = 6370043.69 m 
Radius of ellipsoid at the surface (WGS84) = 6368458.20 m 

6.3 Effect of the GEOID (EGM08) or the angular fluc- 
tuations of the gravitational potential 

The geoid undulation of the geoid reflects the fluctuation in height around the 
surface of the reference ellipsoid (WGS84 here) which defines an equipotential 
surface of gravity. The mean sea level is a good illustration of this geoid 
surface, (cf. 



/Hf it jr Umax H 

GM ^ \ - x - ,a " 



V(r, A,0') = — [1 + ^ ^2 y P n m(sin((j)'))(C nm cos(m\) + S nm sin(m\))} 

(72) 



n=2 m=0 



Normalized associated Legendre function = P nm 
Normalized gravitational coefficients = C nm , S nm 
EGM96 nmax=360, resolution = 1 ° = 60 nautical miles= 111.12 km 
EGM08 spherical harmonic degre and order 2159, resolution= 6/ = 10 nau- 
tical miles= 18.52 km 
Height anomaly = £(r, A, 0') 

7i ir Umax 

(jrlvl \ — -v ^ — -v a n _ _ _ 

C(r, A,0') = [> > (-) P nm (sin((f)'))(C nm cos(mX) + S nm sin(mX))] 

7 (m r z — ' z — 4 r 

1 ^ I n =2 m=0 

(73) 

We have calculated the value of the gravitational potential at the surface 
ellipsoid (WGS84) using the value of the height anomaly in EGM08, thanks 
to the program hsynth_WGS84.f EGM2008 Harmonic Synthesis, provided 
by NGA, Geodesy and Geophysics Basic and Applied Research (|14|) (see 
figure [7]). 

V(r, A, 4>') = -— - C(r, A, 0')t(0) (74) 
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6.4 Bouguer and free air anomaly corrections of the 
gravity inside mountains above the ellipsoid 

There are two additional gravity corrections to take into account: the Free- 
air anomaly correction and the Bouguer anomaly correction. The Free-air 
anomaly is an altitude correction of the gravity. The Bouguer anomaly cor- 
rects the gravity from the presence of mountains that adds some masses close 
to the point of measurement of the gravity. To take this effect into account, 
it is assumed that an infinite horizontal plate of thickness h and density p is 
present between the ellipsoid and the top of the mountain where the gravity 
is supposed to be measured. (cf. [15J & [13J ) . While the measurement is made 
in a mountain above the ellipsoid but below the surface of the mountain, the 
Free-air anomaly and Bouguer anomaly are given as follows: 



Free air gravity anomaly = A gp 
Observed gravity = g 

theoretical gravity at the level of the ellipsoid = 7 

height of the surface of the earth above the ellipsoid = h 

depth of the point under study below the surface of the earth and above the 

ellipsoid = d 



Ag F = g- 1 + AnpGd - (^-)(h - d) - {\){^){h ~ df + Sg A (75) 
AnGpd = 2.238 i(T 6 d (76) 



d'y 
dh 



2-(l + f + m-2f (sin 2 ((f)))) 



(77) 



a 2 b 

m = co 2 -— = 0.00344978650684 
GM 

Angular rotation velocity of earth 



co = 7.292115 10- u rad./s 



Flattening = / 





(78) 
(79) 



dh 2 



Bouguer gravity anomaly = Ag B 



$9b 



2nGph = 1.1195 W~ b h 



(80) 
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Gravitational attraction of Bouguer plates = 5gs 



p = 2.67t/m 3 (81) 
^ = .8710- 5 e- 116 ^ 1000 ) 1 ' 047 ] (82) 

Atmospheric correction = Sg^ 

The height of the Earth surface above the ellipsoid is computed here thanks 
to the GTOPO30 survey which provides all over the world a resolution of 
about half a nautical mile. (Credit: U.S. Geological Survey Department of 
the Interior/USGS) (see figure [S]). 



7 Putting alltogether the corrections (elliptic- 
ity, geoid, Bouguer and Free-Air) in the cal- 
culation of the relativistic potential for the 
CERN to Gran-Sasso neutrino delay time cor- 
rection 

The ellipsoid correction to the General Relativistic delay correction is applied 
in the following way: 

In order to take into account the ellipticity of the earth in the calculation 
of the relativistic gravitational potential $ along the neutrino travel, we 
have calculated for each point along the neutrino travel inside the Earth, the 
value of the ellipsoid radius r e u ipsoid ((f)) and then the distance to the center 
r is normalized to this radius of the ellipsoid at this underground point, 
afterward the relativistic potential is computed to this renormalized point rl. 
Some other method can be used (cf. [16]). 

rl = r —— (83) 

T ellipsoidyr ') 

$(r) $(r/) = T —[1t) ( 84 ) 

r : distance from the center of the Earth to the point under study along the 
neutrino road. 

r eiUpsoid{4') '■ distance from the center of the Earth to the surface of the 
ellipsoid 

k = 1.0000045084305260 : squared ratio between the surface of the ellipsoid 
and the sphere of mean radius. 
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The geoid correction to the gravitational potential is applied in such a way 
that it is maximal at the surface of the ellipsoid and null at the center of the 
Earth, this is a first order modification, of course more refined correction can 
improve this calculation. 

1,^-^ + a <M^ (85, 

For the points of the trajectory of the neutrinos which are above the 
ellipsoid but inside the mountains, we have integrated the gravity acceleration 
with the Bouguer and the Free-Air correction in order to add a correction 
term to the gravitational potential $ used in the Schwartzchild metric. For 
the points located below the ellipsoid we have not applied any correction 
to this relativistic gravitational potential $. The result of the computation 
(Monte Carlo method with 10 9 iterations ) with all these corrections to the 
calculation of the delay effect is St = 4.1863 picosecond which is equivalent 
to a length delay of 1.2550mm, also equivalent to a reduced apparent value 
of the speed of light of ^ = 1.7179 1(T 9 . 



8 Conclusion 

Based on the Schwartzchild metric, a general relativistic metric adapted 
to the Earth interior is detailed here with geophysics density modelization 
(PREM). Non rotation and spherical symmetry is assumed to compute an 
effective zonal prametric general relativistic "gravi-potential" term. The de- 
lay of travel for neutrino crossing the Earth from one point of the surface 
to another distant point is extensively detailed. A Monte-Carlo computa- 
tion of the integral of the general relativistic time delay effect on neutrino 
travel from CERN to Gran-Sasso leads to a St = 4.1888 picosecond in this 
simple spherical irrotational model. This delay effect is also computed and 
fitted for other distances crossing the Earth. If a neutrino crossed the whole 
Earth diameter, a "one inch delay effect" (2.59cm) would be observable of 
St = 86.65 picosecond . Some secondary effects are later included as the el- 
lipticity of the Earth using the WGS84 system , the GEOID gravitational 
potential fluctuations using the EGM08 survey and also the Bouguer and 
Free-Air anomaly gravity corrections using the GTOPO30 survey. The com- 
putation of all these corrections leads to a St = 4. 1863 picosecond delay effect. 
This general relativistic effect cannot account for the St = (57.8±7.8(s£ai.) + 
8.3 — 5.9(sys.))ns seen by the OPERA experiment in Gran-Sasso. Even the 
Icarus experiment St = (0.3 ± 4.0(stat.) ± 9.0s(syst.)ns which does not see 
any discrepancy between the speed of light and the neutrino speed cannot 
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measure such a tiny delay effect due to the general relativistic impact of the 
mass desnity distribution inside the Earth. 
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Figure 1: Density of the Earth interior (PREMC) 
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Gravity acceleration inside earth 
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Figure 2: Gravity acceleration inside the Earth 
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Figure 3: Gravitational potential inside the Earth 
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Figure 4: General relativistic term of gravitational potential — $(r) inside 
the Earth compared to a pointlike gravitational potential 
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Figure 5: Neutrino delay of arrival time versus neutrino baseline length 
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Figure 6: Neutrino length delay of arrival time versus neutrino baseline 
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Figure 7: Height anomaly versus distance of travel from CERN to Gran-Sasso 
Laboratory 
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Elevation above sea level of earth surface along neutrino travel 
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Figure 8: Height of the Earth surface versus distance of travel from CERN 
to Gran-Sasso Laboratory above sea level from GTOPO30 (Credit: U.S. 
Geological Survey Department of the Interior/USGS) 
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